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Knowledge of the dynamical behavior of correlations with no classical counterpart, like 
entanglement, nonlocal correlations and quantum discord, in open quantum systems is 
of primary interest because of the possibility to exploit these correlations for quantum 
information tasks. Here we review some of the most recent results on the dynamics of 
correlations in bipartite systems embedded in non-Markovian environments that, with 
their memory effects, influence in a relevant way the system dynamics and appear to be 
more fundamental than the Markovian ones for practical purposes. Firstly, we review 
the phenomenon of entanglement revivals in a two-qubit system for both independent 
environments and a common environment. We then consider the dynamics of quantum 
discord in non-Markovian dephasing channel and briefly discuss the occurrence of revivals 
of quantum correlations in classical environments. 

Keywords: Open Quantum Systems; Dynamics of Quantum Correlations; Non- 
Markovian Environments. 



1. Introduction 

Realistic quantum systems are important both for their fundamental r olel] and for 
up-to-date technologies^. They are open, in the sense that their correlations with 
surrounding environments are of statistical nature and unavoidable^. The character- 
ization of correlations in composite open quantum systems is essential for quantum 
information and computation tasks^. To a specific property of the system it can 
be made to correspond a type of correlation that is used to define the property it- 
self. In fact, entanglement represents correlations related to non-separability of the 
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state of a composite quantum system; violation of Bell inequalities^^ represents 
quantum nonlocality, that is cor rela tions not reproducible by any local classical 
model; quantum discord (entropicP^ and geometric^) represents quantumness and 
indicators of purely classical correlations can also be defined. Each of these types 
of correlations can be useful for specific tasks. For ex ample , entanglement is the 
main resource for teleportation and su perden se coding^U^, n0 nlocal correlations 
guarantee secure quantum crypto graph y I, entropic discord is a resource for de- 
terministic one-q ubit computatior 13 | 14 | 15 | anc j geometric discord for remote state 



preparatio n * [ The above correlations, each linked to a specific property of bipar- 
tite quantum systems, have been extended to multipartite systems and quantified 
in terms of d istances betwe en the state of the system and its closest states without 
that propert y9 | 18 | 19 | 20 | 21 | ^ - g therefore essential to understand, in a given open 
system, the dynamics of the various correlations and in particular how much and 
for how long they persist. 

The dynamical behavior of correlations present in a composite open quantum 
system strongly depend on the noise produced by the surrounding environment. In 
this sense one of the most important aspects of the environment is if it can be de- 
scribed cither as memoryless (Markovian) or as with memory (non-Markovian) . The 
characterization of the environment as Markovian or non-Markovian, with respect 
to a given system, is determined by the ratio between its typical correlation time and 
the system relaxation time. For two-qubit systems in Markovian environments, some 
kind of correlations may be subject to early-stage disappearance (ESD, or sudden 
death), meaning that they completely disappear at a finite time in spite of exponen- 
tial decay of singlc-qubit coherence. This has been shown, both t heoretically and ex- 
perimentally , in th e case of entanglemenlj 22 l 2 3 l 2 4 l 2 5 | 26 | 27 | 28 | 29 | ^ inequal- 
ity violation o ' '. This phenomenon strongly limits the time when correlations can 
be usefully exploited. Differently, two-qubit systems i n structured non-Markovian 
environment s may pre sent phenomena such as revival js ?2 | 33 | 34 | 35 | 36 | 37 | 38 | 39 | 40 | '| T] 
and trapping32E3E31 of correlations, thus overcoming ESD. The study of the dy- 
namics of correlations under non-Markovian noise thus appears to be fundamental 
for quantum information purposed 2 -^ and it shall be the main subject of this paper. 

In particular, the aim of this review is to summarize some of the most recent 
results on the dynamics of correlations in bipartite qubit systems embedded in 
specific non-Markovian environments. We shall first review the phenomenon of cor- 
relation revivals in a two-qubit system both for independent environments and for 
a common environment. Secondly, we describe how structured environments, as 
photonic crystals, may protect correlations against decoherence. We finally discuss 
the dynamics of quantum correlations in non-dissipative dephasing channels and in 
classical environments given by random external fields. 
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2. Revivals of entanglement 

In this section we review the phenomenon of entanglement revivals for non- 
Markovian bosonic reservoirs, both in the case of independent environments and 
in a common environment. 



2.1. Two-qubit initial states 

We first define the initial states typically chosen for the analysis and describe their 
properties. We consider initially entangled two-qubit states within the class of X 
stated^, whose general structure is, in the standard computational basis B — {|1) = 
|11>,|2> = |10>,|3} = |01>,|4} = |00>}, 



Px 



/ pn p i4 \ 

P22 P23 
P23 P33 
Vpi4 p A4 / 



(1) 



X-struct ured density matri ces may arise in a wide variety of physical 
situationj^|47|48|49|k)|51|52| 

These states are also encountered as eigenstatesin 
all the systems with odd-even symmetry like in the Ising and the XY 

modelsPMl. 

Moreover, for many physi cal dynamics of open quantum systems an initial X struc- 
ture is maintained in time^l This shall occur, in particular, in the following cases 
described in this review. X states contain the class of extended Werner-like (EWL) 
states, defined aM^ 



p* =r|$)($| + (l-r)l/4, p* = r|*)(*| + (l-r)l/4. (2) 
In Eq. (|2J), r is the purity parameter, i the 4x4 identity matrix and 

|*>=o|01) + 6|10), |*) = o|00) + 6|ll), (3) 

are the Bell-like states with a real, b — \b\e n and a 2 + \b\ 2 = 1. For a = ±6 = 



1/V2, Bell-like states become the Bell states |$ ± ) = (|01) ± |10>) /y/2, |* ± ) = 
(| 00) ± 11 1)) A/ 2. The states defined by Eq. ([2]) reduce to the well-known Werner- 
like stateJSSESl 

when their pure part is a Bell state. For r = EWL states become 
totally mixed states, while for r = 1 they reduce to the Bell-like states |$), |\&) of 
Eq. ([3]). One of the important aspects of the EWL states p*(0), /5*(0) is that they 
allow to study the effect, on the dynamics of correlations, of both the initial state 
mixedness and the degree of entanglement of their pure part. 

The ent ang lement of the bipartite system shall be hereafter quantified by the 
concurrence^. For a X state, as defined in Eq. (JlJ , concurrence can be easily com- 
puted and is given bjS^l 

C Px = 2max{0, K\, K 2 }, K\ = |p 2 3| - y/PnPu, K 2 = \pu\ - VP22P33- (4) 

EWL states of Eq. © give equal value for the concurrences, given by 

Cp = Cf = 2max{0, (o|6| + l/4)r - 1/4}. (5) 
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Fig. 1. Schematic of the two-qubit systems with independent environments (left figure) and 
common environment (right figure). 



From previous equation one finds that there is entanglement when r > r* = 
(1 + 4a|fe|)- 1 (for a = \b\ = l/y/2 one has r* = 1/3). In particular, the Bell- 
like pure states |$),|^) have the same degree of entanglement (concurrence) 
C$(0) = C*(0) = 2as/l - a 2 . 



2.2. Independent bosonic environments 

We consider a system composed by two independent parts each made of a qubit 
placed inside a zero-temperature bosonic environment (see left part of Fig. Q}. The 
two qubits A and B, with ground and excited states |0), |1) and transition fre- 
quency ujq-, are identical and initially entangled. A procedure has been introduced, 
applicable to a general multipartite system made by independents subsystems, that 
allows to obtain the reduced density matr ix of t he system at the time t by the 
knowledge of that of the single subsystems^ESl. Here we summarize this proce- 
dure for the case of a two-qubit system. For independent parts the total time evo- 
lution operator can be factorized in two terms relative to the two parts. Under 
these conditions, given the dynamics of each qubit as p^,(t) = J2 W A^,(t)p^,,(0), 
Pfj' (*) = Emm' B jj' n ' (*)/°mm'(0), the two-qubit density matrix is simply given byS2J 

PwMt)= E A%{t)B]p\t)p w , mm ,{Q)> (6) 

W ,mm f 

where the indexes i, j, l,m = 0, 1. This procedure is ea sily generalizable to an arbi- 
trary multipartite system of qudits (d- level subsystems)^ and it can be also applied 
to the case of two qubits in a common environment provided that the qubits are 
sufficiently separated in order to neglect the dipole-dipole interaction and their 
distance is larger than the spatial correlation length of the reservoir. 

The dynamics of each qubit coupled to the environment modes is described by 
the spin-boson HamiltoniarP 

H = hujQ(j + <T- + t^kb\b k + ^ HdhV+bk + g* k (J-b\), (7) 

k k 

where ojq is the qubit transition frequency, <j± are the qubit raising and lowering 
operators, the index k labels the environment field modes with frequencies u>k, b\ 
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Dynamics of quantum correlations in two-qubit systems within non-Markovian environments 5 

and bk are the creation and annihilation operators, respectively, and gu the coupling 
constants. In the following, we shall give general expressions for the concurrence 
valid for any local qubit-environment interaction described by the Hamiltonian of 
Eq. (JT]) . When the environment is at zero-temperature and the qubit is initially in a 
general superposition state of its two levels, the singlc-qubit reduced density matrix 
p s (t) (S = A, B) takes the forirP 

/pfi(o)|g(*)l 2 pfo(ok(t) \ 

P S (t) = . (8) 

V P$i(0)Q*(t) Po S o(0)+pfi(0)(l-kW| 2 )/ 

The singlc-qubit dynamics thus depends only on the function q(t) (0 < \q(t)\ < 1, 
q(0) = 1), whose explicit time-dependence contains the information on the envi- 
ronment spectral density and the coupling constants. Using Eq. ([6]). one can then 
determine the two-qubit dynamics. In particular, in the standard basis B = {|1) = 
|11), 1 2) ee 1 10), |3) ee |01), 1 4) ee |00)}, the diagonal elements of the reduced density 
matrix p(t) at time t result to be 

= Pn(0)k(t)| 4 , pa® = Pll (0)k(t)| 2 (l - \q(t)\ 2 ) + P 22(0)W)\ 2 : 
pas(t) = pu(0)| g (t)| a (l ~ k(t)\ 2 ) + p 33 (0)k(i)| 2 , 

P4A® = 1 - [pii(t) + p 22 (t) + p 33 {t)}, (9) 
and the non-diagonal elements 

Pu(t) = p 12 (0)q(t)\q(t)\ 2 , p 13 (t) = p 13 (0)q(t)\q(t)\ 2 , p u (t) = p u (0)q(t) 2 , 

p 23 (t) = P 2 3 (0)\q(t)\ 2 , P2 4 (t) = Pis(0)g(t)(l - \q(t)\ 2 ) + PM(0)g(t), 

P3i(t) = p la (0) g (t)(l - \q(t)\ 2 ) + pu(0)q(t), (10) 

with pij{t) = Pji(t), p(t) being a hermitian matrix. Eqs. ([9]) and (fT0|) give the 
two-qubit density matrix evolution for any initial state dependent only on q(t). 

This way, for the case of the above spin-boson Hamiltonian the concurrence, or 
any other measure of correlations calculated on the system state, can be shown to 
depend only on the initial state conditions and on the single-qubit parameter q(t), 
whose explicit form is determined by the specific spectral density of the environment. 
This explicit dependence of the concurrence on q(t) has been already reported in 
the case of initial Bell-like states of Eq. (J3j> , when the system state remains of this 
form during the evolution (as happens in the case here considered). In particular, 
these expressions are 

C7#(t) = 2max{0,a|&||< ? (£)| 2 }, C*(t) = 2max{0, \b\\q(t)\ 2 [a - \b\(l - \q(t)\ 2 )}}. (11) 

The time-dependent forms of concurrences above show that, if the two-qubit system 
starts from the Bell-like state |$), ESD can never occur being C$(t) proportional to 
a power of the single-qubit coherence |<?(£)|. For initial Bell (maximally entangled) 
states, |'I >± ) and |'I' ± ), the concurrences at the time t take the simple form C$± (£) = 
|<7(£)| 2 , Cy±(t) = |g(i)| 4 , from which it is readily seen that for both these initial 
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states ESD never occurs. Now we shall extend the expressions of concurrence in 
terms of q(t) also for the rather more general class of initial EWL states of Eq. ([2]), 
whose Bell-like states are a subclass. 

During this dynamics, the X-structurc of the initial density matrix is maintained 
and this is the case also for EWL states. The explicit expressions of concurrences 
at time t are given for the two initial states p* and of Eq. ([2]), respectively, by 



C p *(t) = 2max|0,| g (t)| 5 
C p *W = 2max(o,|^)| 2 



ra|6|-iy(l-r) (l - \q(t)\* + \q{t)\^— 
™\b\ - ^-(\q{t)\ 2 {l + 3r - ira 2 ) + Ara 2 - 2(1 + r)) 



(12) 

where \b\ = yl — a 2 . These equations are quite general since their form does not 
depend on the particular choice of the environment, but only on the Hamiltonian 
model of Eq. ([7]) and on the chosen initial state. Once the environment structure is 
specified, the explicit form of q(t) is obtained that in turns determines the explicit 
time-dependence of concurrence. 

2.2.1. Bosonic environment with Lorentzian spectral density 

The Hamiltonian of Eq. ([?])■ which may describe various spin-boson systems, is now 
taken to represent a qubit formed by the excited and ground electronic state of 
a two-level atom interacting with the quantized modes of a high-Q cavity. This 
dissipative qubit- environment interaction gives the so-called amplitude dam ping 
channel. At zero temperature, this open quantum system has an exact solution^? 
In the case of a single excitation in the atom-cavity system, the effective spectral 
density J(w) is taken as the spectral distribution of an electromagnetic field inside 
an imperfect cavity supporting the mode u>o, resulting from the combination of the 
reservoir spectrum and the system-reservoir coupling, with F related to the micro- 
scopic system-reservoir coupling constant. This spectral density has the Lorentzian 
forrrP 

JM = (rA 2 /27r)/[( Wo -^ 2 + A 2 ], (13) 

where A, defining the spectral width of the coupling, is connected to the reservoir 
correlation time tb by the relation tb ~ A -1 . In fact, it can be shown that the 
reservoir correlation function, corresponding to the spectral density of Eq. (|13[) , 
has an exponential form with A as decay rate. On the other hand the parameter 
r can be shown to be related to the decay of the excited state of the atom in the 
Markovian limit of flat spectrum. The relaxation time scale tr over which the state 
of the system changes is then related to F by tr « T -1 . The exact solution for the 
single-qubit density matrix (at zero te mpera ture) with the spectral density given 
by Eq. flT5J} has the form of Eq. © wittP^ 

\q(t)\ = c - xt/2 [cos((ft/2) + (X/d) sin(dt/2)] , (14) 
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Fig. 2. Concurrence as a function of Ft for initial Bell-like states |<t>) (blue dashed line) and 
(orange solid line) for the case of independent non-Markovian environments with A/T = 0.01. The 
red solid curve on the left represents the concurrence for | VP) under Markovian condition A/T = 3. 
The initial Bell-like states have a 2 = 1/3, 



where d = \/2TX — A 2 . A weak and a strong coupling regime can then be distin- 
guished for the single qubit dynamics. The weak regime corresponds to the case 
A/T > 2, that is tr > 2tb- In this regime the relaxation time is greater than the 
reservoir correlation time and the dynamics is essentially Markovian with an expo- 
nential decay controlled by T. In the strong coupling regime, that is for A/T <C 2, 
or tr <C 2tb , the reservoir correlation time is greater or of the same order of the re- 
laxation time and non-Markovian effects become relevant. Substituting the explicit 
expression of \q(t)\ of Eq. (TT4"|) in Eqs. © and (|10[) . the two-qubit density matrix 
at time t and its concurrence are obtained for both the initial states of Eq. 

The time behavior of the concurrences C^,(t) and C*(i) as a function of the 
dimensionless time Tt are plotted in Fig . [2] . For Markovian conditions (A/T > 2) 
the phenomenon of ESD is retrieved^, that is entanglem ent disappears at a fi- 
nite time, as also observed both in an all-optical experimenlPSl and in two atomic 
ensembles^!. When non-Markovian effects are relevant (A/T <^ 2), revivals of entan- 
glement may appear after finite periods of time when disentanglement is complete. 
This condition of non-Markovianity can be achieved, for example, in cavity QED 
experimental configurations by using Rydbcrg atoms with lifetimes T at ~ 30ms, 
inside Fabry-Perot cavities with quality factors Q ~ 4.2 x 10 10 corresponding to 
cavity lifetimes^' T cav w 130ms: these values indeed correspond to 2A/T ~ 0.2. The 
phenomenon of revivals is induced by the memory of the reservoir, which allows the 
two-qubit entanglement to reappear after a dark period of time, during which the 
concurrence is zero. These revivals, although they effectively extend the possible 
usage time of entanglement, decrease with time and eventually disappear after a 
certain critical time. In conclusion, revivals of entanglement, after finite periods of 
entanglement disappearance, are linked to the qubit-environment non-Markovian 
dynamics and, after their discovery, they have been explained in terms of transfer 
of correlations back and forth from the two-qubit system to the various parts of 
the total system. This correlation transfer can be seen as due to the back-action 
via the environment on the system, which creates correlations between qubits and 
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environments and between the environments themselves: in par ticular, co rrelations 
may build up between the two independent quantum reser voir d ^ 2 1 6 3] (^his phc- 
nomenon has been named sudden birth of entanglement in reservoirs). In this sense, 
the nature of the entan glemen t revivals described above, occurring for bosonic in- 
dependent environmentally, differs from that of revivals obtai ned for interacting 
qubits embedded in a common Markovian bosonic reservoiJ^GHS where revivals are 
just due to a direct dipolc-dipole interaction among. Collapse and revival features of 
the entanglement dynamics of different polarization-entangled photon states have 
been observed in an all-optical experimental setup, when one of the travelling pho- 
tons passes through a non-Markovian quant um e nvironment simulated by a suitable 
Fabry-Perot cavity followed by quartz plates^l. It has been also recently seen that 
only ESD with no revivals may occur in the case of independe nt non -Markovian 
baths with broadband spectra typical of solid-state nanodevice d^®^ ! Recently it 
has been shown that the presence of the system-bath coherence in the initial state 
may change dramatically the entanglement dyna mics of two qubits in the pres- 



ence of a non-Markovian quantum-mechanical battpSl. Entanglement dynamics of 
a system of N qubits interacting with N independent reservoirs in both Markovian 
and non-Markovian regimes has been also investigated. Entanglement revivals of 
qubits at instantaneous points of disappearance or after a finite interval of abrupt 
va nish ing due to the memory effect of non-Markovian reservoirs have been pointed 

Also the d ynamic s of nonlocal correlations, identified by violations of the CHSH- 
Bell ine 

quality™ 

has been studied in independent non-Markovian bosonic en- 
vironments, like the ones described here, and compared with the dynamics of 
entanglemenljSSI. Revivals of Bell inequality violations have been thus found but 
they show up for quite strong non-Markovian features. In particular, there may ex- 
ist time regions when the Bell inequality is not violated even in correspondence with 
high values of concurrence (C ~ 0.8: this indicates that even highly entangled sys- 
tems cannot be exploited with certainty in contexts where the presen ce of n onlocal 
correlations is required (for instance, in secure quantum cryptogra phyinni). 

An all-optical experimental setup has been recently employed to observe and 
engineer entanglement oscillations of a pair of polarization qubits in an effective 
non-Markovian channel. Programmable oscillations of entanglement are achieved, 
where the entangled state obtained at the maximum of the revival after sudden 
death violates Bell inequality by 17 standard deviations^!. 



2.3. Common bosonic environment 

We now review the main features of entanglement dynamics when two nonintcr- 
acting qubits are embedded in a common non-Markovian bosonic environment (see 
right part of Fig. 1). As the literature of this to pic is extensive, we will focus here 
on a simple model of dissipative dynamic d^^ . 

The Hamiltonian for the total system, in the dipole and the rotating-wave ap- 
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proximations, and in units of h = 1, can be written as 

H = uxa^a^ + uj 2 af a {T) + ]T u k b\b k + (a^ + a 2 (jf ) ) £ 9kb k + h.c.(15) 

fc fc 

where b k , b k are the creation and annihilation operators of quanta of the reservoir, 
a± and uij are the inversion operators and the transition frequency of the j-th 
qubit (j = 1,2), respectively, uj k is the frequency of the reservoir fc-th mode, and 
ctjg k describe the coupling strength between the j-th qubit and the fc-th mode of 
the reservoir. 

Here, cxj are dimensionless real coupling constants measuring the interaction 
strength of each single qubit with the reservoir. In the case of two atoms inside a 
cavity, e.g., different values of ctj can be obtained by changing the relative position 
of the atoms in the cavity field standing wave. We denote with cut = (a? the 
collective coupling constant and with r\- = otj/ocx the relative interaction strength. 



2.3.1. One- excitation dynamics 



We begin by considering the case in which only one excitation is present in the 
system and the reservoir is in the vacuum. We consider the initial state 



*(0))= Coi|l)l|0> 2 + Co2[0)l|l>2 



10 



k R, 



(16) 



where coi and Co 2 are complex numbers, \0)j and (j — 1, 2) are the ground and 
excited state of the j-th qubit, respectively, and \0 k ) r is the state of the reservoir 
with zero excitations in the fc-th mode. 

The time evolution of the total system is given by 



|*(t)> - ci(t)|l)i|0) a |0) fl + c 2 (t)|0) 1 |l) 2 |0) fl + X)cfc(*)|0)i|0) 5 



k/R, 



(17) 



|lfc)ii being the state of the reservoir with only one excitation in the fc-th mode. 

In the standard basis, the reduced density matrix, obtained from the density 
operator |^ f (f))(*(f)| after tracing over the reservoir degrees of freedom, takes the 
form 



P(t) 



/0 

| Cl (i)| 2 Cl (t)c* 2 (t) 
0cl(t)c 2 (t) |c 2 (t)| 2 

\o 



i-N 



(18) 



The two-qubit dynamics is therefore completely characterized by the amplitudes 
Ci j2 (i)- For certain specific structures of the reservoir one can obtain the exact 
analytical expressions of ci j2 (t) by the Laplace transform method. 

We consider a structured reservoir describing the electromagnetic field inside 
a lossy cavity. This case can be modeled by a Loren tzia n broadening of the fun- 
damental mode cavity, as given by Eq. ([T3|) . In RefP^ the exact, and therefore 
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non-Markovian, analytical expression for the amplitudes C\^{t) has been derived. 
Using these quantities, one can readily calculate the time evolution of entan- 
glement, as measured by concurrence, which for our system is simply given by 
C(t) = 2 |ci (i)c2 (t) | ■ An interesting feature of the dynamics is the existence of a 
non-zero stationary value of C due to the entanglement of the decoherence-free state 
1^-}. The maximum stationary entanglement C™ ax ~ 0.65 is obtained for initially 
factorized states. 

For weak couplings and/or bad cavity and for an initially separable state, the 
concurrence increases monotonically up to its stationary value; whereas, for initially 
entangled states, the concurrence first goes to zero before increasing towards C s . 
The strong coupling/good cavity case is more rich and presents entanglement oscil- 
lations and revival phenomena for every initial atomic states. One can prove that 
for maximally entangled initial states the revivals have maximum amplitude when 
only one of the two atoms is effectively coupled to the cavity field, i.e. for n = 0, 1. 
This revivals are similar to those present in the independent environments case dis- 
cussed in Sec. 12. 2\ as they are also due to the reservoir memory effect. However, 
in the case of common environment the revivals are enhanced and the amount of 
revived entanglement is huge, being comparable to the previous maximum. 

Based on this dynamics, a detailed analysis of the experimental implementation 
of cavity-mediated entanglement generation in the context of trapped ion cavity 
QED has been presented in RefM-^ 1 . The rich dynamics of the system, showing 
quantum beats and other tru ly nonclassical phenomena, has been studied in the 
off-resonant regime in RefP^l. Finally, it has been demonstrated that the loss of 
entanglement can be completely suppressed by a series of app ropriate projective 
measurements, a consequence of the quantum Zeno 



eries otappr 
efifecPaZal. 



2.3.2. Two-excitation dynamics 

Let us now consider the case in which two excitations arc initially present in the sys- 
tem. In this case the derivation of the exact analytical solution is more complicated. 
For the sake of simplicity, we assume that two identical atoms inte ract resonantly 
with a Lorentzian structured reservoir with equal coupling. In RefP it has been 
shown that, by using the pseudomode approach^! an d establishing a connection 
with a three-level ladder system, it is still possible to solve the dynamics without 
performing any approximations. Armed with the exact solution we can prove that 
ESD and entanglement su dden b irth (ESB) phenomena occur in the non-Markovian 
common reservoir scenarid ^^H 

Entanglement sudden death and entanglement revivals are basically due to two 
combined and intertwined effects: the backaction of the structured reservoir and the 
reservoir-mediated interaction between the qubits. In order to understand the role 
played by each of these effects it is useful to compare the dynamics with the cases of 
the common Markovian reservoir, and two independent non-Markovian reservoirs. 

In a common Markovian reservoir a period of complete disentanglement is fol- 
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lowed by a revival of entanglement, but no oscillations are present. Such a revival 
is due to the action of the common reservoir which tends to create quantum cor- 
relations between the qubits, providing an effective coupling between them. The 
feedback of information from the reservoir into the system, characterizing the non- 
Markovian dynamics, enhances the appearance of ESD regions, since it tends to 
recreate the conditions that led to the first ESD period. 

A comparison between the non-Markovian independent and common reservoir 
cases reveals that, for the same type of reservoir spectrum, ESD (dark periods) 
regions are much wider in the independent reservoirs case than in the common 
reservoir case. Since both cases take into account memory effects, this suggests that 
the reservoir-mediated interaction between the qubits, in the common reservoir 
scenario, effectively counters the fast disappearance of entanglement. 

Non-Markovian effects also influence strongly the dynamics for initially factor- 
ized states, when the reservoir-mediated interaction between the qubits leads to 
entanglement generation. While, for independent reservoirs, sudden birth of entan- 
glement does not appear, for the common reservoir case, ESB does occur. Moreover, 
in a structured reservoir such a phenomenon presents new interesting features, com- 
pared to the Markovian case. In general, the reservoir memory prolongs the initial 
disentanglement. In this case, therefore, the reservoir backaction dominates over 
the tendency of the common reservoir to create entanglement between the qubits. 
ESB periods and disentanglement revivals become more frequent and numerous for 
stronger non-Markovian conditions. 



3. Entanglement trapping in independent structured environments 

Strategies for preventing the decay of quantum correlations, initially present in a 
composite quantum system, and thus to fight decohcrence are important for effec- 
tive realization of quantum computers. Some of these strategics, applied to a single 
qubit, are quantum error correction^, dynamical deco upl ing of the quantum sys- 
tem from its environment^^] or quantum Zeno effedP^. All of these procedures 
require external controls. It has been then shown that ent angl ement protection is 

loni 

possible by dynamical decoupling base d on n onlocal pulses 2 ^ and by Zeno effect 



on qubits in a common bosonic reservoi J^^^^ i On the other hand, another possible 
way to reduce the detrimental effect of the environment, without using external 
controls, could consist in embedding the qubits in suitably structured environments 
where, for example, population trapping is achievable. Th is is k nown to happen 
in photonic band-gap (PBG) materials or photonic crystalJSHSS. Here we briefly 
review the re sults obtained for two noninteracting qubits in independent photonic 
crysta lplll. 

Each environment is assumed to be a zero-temperature three-dimensional pe- 
riodic dielectric with isotropic photon dispersion relation Uk (photonic crystal)^. 
The Hamiltonian model is always of the spin-boson type as in Eq. (0 with a noise 
spectral density typical of a PBG material. By symmetrizing w^, one produces pho- 
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Fig. 3. Asymptotic values of concurrences C$(t = oo) (orange solid line) and C\n(t = oo) (green 
dashed line) as a function of (5//9, starting respectively from the initial Bell state l"!^} and l^^). 

tonic band gaps at the spheres |k| = mir/L (m = 1,2,.. .), where L is the lattice 
constant. In such ideal photonic crystals, a PBG is the frequency range over which 
the local density of electromagnetic states and the decay rate of the atomic pop- 
ulation of the excited state vanish. Near the band-gap edges the density of states 
becomes singulaJ^, the atom-field interaction becomes strong and one can expect 
modifications to the spontaneous emission decay. For k = tt/L the dispersion rela- 
tion near the band edge frequency ui c can be approximated by Wfe = oj c + D(k — k) 2 , 
where^l D = to c /k 2 . This dispersion relation is isotropic since it only depends on the 
magnitude k of the wave vector k. In the case of a two-level atom, with atomic tran- 
sition frequency o>o, embedded in such a m ater ial, the explicit form of singlc-qubit 
coherence q(t) = f(5, (3, t) has been obtained^! in terms of the detuning 5 = uj —uj c 
and of a parameter (3 defined as {3 Z I 2 = uj^ 2 d 2 /6Treohc 3 , where eo is the Coulomb 
constant and d the atomic dipole moment. The parameter /? 3 / 2 , being proportional 
to WqG? 2 , represents the strength of the atom-reservoir coupling. The time evolution 
is scaled by (3 and for large times ((3t ^> 1) one has |g(£)| 2 — ► constant, that is a 
population trapping 2 ^. 




The dynamics of entanglement is obtained by the expressions of the concurrences 
in terms of q(t) = f(S, (3, t), as for example those of Eq. (fTT|) when the two qubits are 
prepared in a Bell-like state. The concurrence is nearer to its maximum value 1 for 
atomic frequencies farther from the band edge and deeper inside the gap. In Fig. [3] 
the asymptotic values of C$± (t —¥ oo) and Cj± (t — > oo) are plotted as a function 
of the ratio S//3 and it can be appreciated that C$±(t —J- oo) is slightly larger than 
C^±(i — > oo). We see that, when 8//3 < 0, these asymptotic values differ from zero 
and decrease very rapidly going near the edge. It is evident that the more the atomic 
transition frequency is far from the band edge and inside the gap (S//3 < 0), the 
higher is the preserved entanglement. An isotropic disp ersion relation is known to 
lead also to the formation of atom-photon bound stateiP^EiS that in turn lead to 
entanglement trapping near the edge even for uiq outside the band gap (6/(3 > 0), 
although with small asymptotic values (< 0.4; sec Fig. [3J. For an initial Bell-like 
state one finds that smaller initial entanglement leads to lower asymptotic value of 
concurrence and that the entanglement trapping phenomenon does not depend on 
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the initial value of a. It is worth to point out that, in a zero-tcmperaturc Markovian 
environment, th e initial Bcll-likc state |\&) of Eq. (|3|) would undergo entanglement 
sudden deatfPSlfor a < 1/V2; differently, in a PBG material entanglement sudden 
death and its consequent limit on the entanglement usage can be prevented^!. 

The results described above arc valid for ideal PBG materials, thus one may 
ask to what extent they hold for real crystals. In fact, in real crystals with finite 
dimensions a pseudogap is typically obtained with a density of states, although not 
exactly zero, strongly smaller than that of free space. For quantum dots inside a 
real PBG material, for example, a pronounced inhibition of spontaneous emission 
up to ~ 30% has been experimentally observed^ which would give, on the basis 
of the relation between population and concurrence given by Eqs. (fT2"|) and (jTTJ) , a 
partial inhibition of entanglement decay. 

Other analyses have been done, after the above one^^l, f or two-qubit entangle- 
ment in one-dimensional p hotoni c crystals^! in a common PBG material and with 
a dipole-dipole interactio n 1 ^^ , where sudden death and sudden birth of entan- 
glement have been found. We also point out that nonlocality, identified by Bell 
inequality violations, can be also efficiently protected in photonic crystals^. The 
dynamics of quantum discord has been also recently studied in a common PBG 
material, finding that the quantum discord can maintain a constant value in the 
long-time limit even when entanglement suddenly disappeari"" 



4. Dynamics of correlations other than entanglement 

Here we consider the dynamics of quantum correlations, identified by quantum dis- 
cord, under different conditions. Relations among quantifiers of entanglement, non- 

IQQI 

locality and purity have been studied in dynamical contexts^ 1 to better understand 
if and how these quantum properties can be connected to each other. Moreover, it 
has been shown that bo th q uantum and classical correlations can be protected via 
the quantum Zeno effect!^. 



4.1. Frozen discord in non-Markovian dephasing channels 

In the last years a huge deal of attention has been devoted to the dynamics of 



correlatio ns other tha n entanglement in presence of both MarkoviarPS and non- 
MarkoviarP14014119U 

quan tum environments. For a comprehensive review on quan- 
tum discord we refer to RefP. In general quantum discord is more robust than en- 
tanglement against decoherence, and it does not present sudden death phenomena. 

A peculiar aspect of the dynam ics of quantum discord is that it may remain 
frozen for finite time interval d^' ^^ eve n in presence of decoherence, as also proven in 
an all-optical experimental setup^^l In fact, it has been very recently theoretically 
demonstrated that, for non-Markovian dephasing, quantum discord can be eternally 



frozen or time invariant)^. This new phenomenon is indeed quite unique, as discord 
is the only nonclassical property which may remain constant in time in presence 
of decoherence even when the state of the system evolves due to the interaction 
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with the environment. If, as several studies strongly indicates, quantum discord is a 
resource for certain quantum technologies, its time invariant properties might play 
an important role. 

The physical system under investigation consists of two qubits under local de- 
phasing channels. Let us assume that the qubits, A and B, are prepared in a state 
of the form 

PAB = \(lAB + Y^c i< ^afj, (19) 

where is the Pauli operator in direction j acting on A(B), Ci is a real number, 

such that < \a\ < 1 for every i, and Iab is the identity operator of the total 
system. For a dephasing channel subjected to either white or colored noise the 
form of the state is maintained during the evolution. In the case of white noise 
the Markovian dissipator takes the form C[p^B)] — l^f PA(B) u f ~ Pa(B)]I^^ 
with 7 the phase damping rate, and j = 1, 2, 3 for the bit, bit-phase, and phase flip 
cases, respectively. For certain classes of states the discord does not decay at all, until 
a finite transition time t. The system indeed experiences two different dynamical 
regimes: at the beginning the state of the qubits evolves in a way to lose only classical 
correlations while quantum correlations are not touched in any way by decoherence. 
This is what we call the classical decoherence regime. In the second phase, after a 
sudden change in quantum and classical correlations, quantum correlations are lost 
while classical correlations are preserved. This is the quantum decoherence phase. 
This effect was named sudden transition from classical to quantum decoherence. 

When the two quits are subjected to local colored noise dephasing channels one 
needs to use a memory- kernel master equation to study the dynamics of each qubit 
subjected to random telegraphic noised In this case, for certain states and regimes 
of parameters, one can distinguish two different behaviours for the dynamics of the 
quantum correlation: the frozen behaviour, and the sudden change dynamics. In 
the first case the dynamics displays multiple sudden transitions: as in the Marko- 
vian case, classical correlation at first decreases while the discord remains constant 
then at a sudden transition point classical correlation becomes constant and dis- 
cord oscillates, and then again discord becomes constant while classical correlation 
oscillates. In the second case first both classical and quantum correlations decrease 
until, at a sudden change point, classical correlations become constant while quan- 
tum correlations exhibit a discontinuous change in the amplitude of the damped 
oscillations. At the following sudden change point classical correlation starts again 
to oscillate and quantum discord changes back to its previous rate of oscillation. 
This behavior reflects the presence of memory effects characterizing the dynamics. 

4.2. Revivals of quantum correlations without back-action 

Here we treat the case when environments arc classical, where back-action is absent. 
One would expect that if the quantum environment of a quantum system is modelled 
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as classical, for example a classical light field coupled to a quantized atom, then 
there might be qualities of the dynamics that this model is unable to describe. 
Intuitively, one might expect that since a classical environment should not be able 
to store quantum correlations in the same way as a quantum environment, the 
capability to describe revivals of quantum correlations might be affected. Similarly, 
if the environment has no back-action on the system, then it would seem that this 
could also affect the correlation dynamics. Ne verthel ess, revivals of entanglement 
can occur for random classical telegraph nois d^^ , for atoms subj ect cither to 
local phase-noisy laserPl or to random external classical fields™ 

It is clearly 

important to exactly know how correlation dynamics behaves when an environment 
is modelled as classical instead of quantum, o r whe n back-action on the system is not 
present. We now briefly review a recent resullP^ that tries to explain how and why 
revivals of correlations, including quantum correlations, can generally occur also if 
the environment is classical, and when no back-action is present. We illustrate our 
discussion with a model characterized by the absence of correlations induced by 
back-action, where the element of randomness is introduced in a very simple way. 

The example we consider is a pair of independent qubits driven by single 
classical field modes with random phase, where it has been shown that revivals 
of correlations do occur both for enta nglement, for the so-called quantum dis- 
cord, and for classical correlation J ^ The dyna mical m ap for the single qubit 
S = A, B is of the random external fields and can be written as 



A s (t, 0)ps(0) = \ X)- =1 Uf (t)ps(0)Z/f f (i), where Uf(t) = e ~ iH * t / h is the time evo- 



lution operator with Hi = i%(o"_|_e _1 *' < — cr-e 1 ^) and the factor 1/2 arises from the 
equal field phase probabilities of the model (more in general, there is a probability 
pf associated to Uf). Each Hamiltonian Hi is expressed in the rotating frame at the 
qubit- field resonant frequency u). In the basis {|1), |0)}, the time evolution operators 
Uf(t) have the matrix form 



where i = 1,2 with qbi = and 02 = tt- This simple model is depicted in the left 
part of Fig. HJ where the two local environments are realized by two single classical 
field modes of the same frequency and amplitude, but each passing across a random 
dcphaser such that the phase of the field at the location of each qubit is either 
zero or n, with probability p = 1/2. The interaction between each qubit and its 
local field mode is assumed to be strong enough so that, for sufficiently long times, 
the dissipation effects of the vacuum radiation modes on the qubit dynamics can be 
neglected. The dynamics is cyclic, as follows from the unitary matrix of Eq. (|20j) (for 
each i = 1,2) becoming the identity for times t = kir/g (k = 0, 1,2, . . .). The same 
initial state is thus retrieved at these times. If one takes into account dissipation 
effects on the system, the dynamics would not be cyclic anymore, but it is expected 
that for dissipation weak enough the qualitative behavior of the dynamics would 
remain the same until a characteristic time. 




(20) 
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Fig. 4. Left figure. Sketch of the physical system. The two qubits are initially entangled. Each 
random dcphascr (RD) is such that the field phase at the location of each qubit is either zero 
or 7r with equal probability, p = 1/2, corresponding to either EjJ or (S = A,B). Right 
figure. Dynamics of entanglement E (orange solid line) and discord D (blue dashed line) for an 
initial Bell-diagonal state Pb(0) = 0.9|<t> + }(<t> + | + 0.1|<f>~)(«l> _ |. D vanishes only at given times 
t n = (2n — l)7r/(4<;) (n = 1, 2, . . .), while E disappears for finite time intervals. Figures reproduced 
from [R. Lo Franco et al, Phys. Rev. A 85, 032318 (2012)]. 



The global dynamical map A applied to an initial state p(0) of the bipartite 
system, p(t) = A(t,0)p(0), is again of the random external fields type, that is, 

<>(*) = I E Uf{t)Uf{t)p{Q)Ut*Wf\t). (21) 

This map A moves inside the class of Bell-diagonal states, that are mixture of 
the four Bell states. Given a state p, we adopt the definitions^^ of entanglement 
E(p) = S(p\\a p ) and quantum discord D(p) = S(p\\x P ) = S(x P ) — S(p), where 
S(p\\a) = — Tr(plogcr) — S(p) is the relative entropy, S(p) = — Tr(plogp) is the von 
Neumann entropy, Xp an d c p are, respectively, the classical state and the separable 
state closest to p. These quantities can be analytically calculated for Bell-diagonal 
states. Right plot of Fig. [4] displays the dynamics of entanglement E and discord D 
starting from the Bell-diagonal state pb(0) = 0.9|$ + ) ($ + | + 0. 1 1 <E> ) ($~ |. Collapses 
and revivals of entanglement occur, as is known to happen for indepe nde nt non- 
Markovian bosonic environments acting as amplitude damping channels^! Quan- 
tum discord also decreases and increases during the dynamics, vanishing only at 
given times. In presence of small dissipation, the same qualitative behavior of cor- 
relations is expected, that is, partial revivals should still occur when the dissipation 
effects are weak enough. 

The above revivals of quantum correlations, as characterized by entanglement 
and discord, occur without transfer of correlations induced by back-action from sys- 
tem to environment. The absence of back-action implies that the initial two-qubit 
correlations cannot be "temporarily stored" as correlations entirely within the en- 
vironments. This differs from the point of view where revivals are interpreted in 
terms of transfer of correlations induced by back-action which correlate the inde- 

pende nt non-Mar k ovian quantum environments, which in turn again re-correlate the 
qubit j32|6l|62|63| 

In our example this mechanism is absent. These revivals can be in- 
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stead explained i n term s of correlations in a classical-quantum state between qubits 
and environments^^. The environments itself play an important role through the 
classical record they keep memory of what operation has been applied to the qubits. 
Finally, it has been also shown that these revivals can be interpreted by referring 
only to the intrinsic features of the two-qubit dynamics, by sho wing that they are 
connected with the increase of a quantifier of non-Markovianit y ' of the singlc- 
qubit dynamics. 

5. Conclusions 

The reduced dynamics of a system interacting with a non-Markovian environment 
can be characterized by relevant memory effects. In this non-Markovian regime, 
correlations present in a composite system may evolve in ways which can be ex- 
ploited in applications requiring their presence. In this paper, we have reviewed 
some recent results concerning the characterization of different kinds of correlations 
in non-Markovian environments. 

In a first part of this review, we have discussed the phenomenon of entanglement 
revivals. In the case of two non-interacting qubits initially prepared in an entangled 
state and subject to local Markovian noise, in front of an exponential decay of 
single qubit coherences, entanglement may disappear in a finite time. Differently, 
in the case of non-Markovian environments, it has been shown theoretically and 
experimentally that entanglement may revive also after finite time intervals when 
it is absent, because of the memory effects in the two-qubit dynamics. 

Following this result, as a possible strategy for preventing the decay of quan- 
tum correlations, we have discussed on the possibility to exploit at the maximum 
non-Markovian effects to protect entanglement from noise. In the case of qubits 
embedded in photonic crystals, entanglement may be preserved for arbitrary long 
times, phenomenon due to the direct link between the dynamics of entanglement 
and that of the excited state population of each qubit; indeed in a photonic-band 
gap material, spontaneous emission of a qubit may be inhibited and population 
trapping may occur, as observed for example in the case of quantum dots. 

In the final part we have also discussed recent results concerning the occurrence 
of revival of correlations also in the case of classical environments such that back- 
action on the system is absent. The mechanism at the basis of these revivals appears 
to be connected to the memory effects in the reduced dynamics and is linked to a 
subtle mechanism that the environment itself can play through the classical record 
they keep memory of what operation has been applied to the qubits. This occurs 
even if the classical environments cannot store any quantum correlations on their 
own, and they do not become entangled with their respective quantum systems. 
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